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I .  INTRODUCTION 

The  problem  of  electromagnetic  coupling  to  a  wire  through  a  small 
aperture  in  an  infinite  conducting  plane  has  been  studied  in  the  past. 
Kajfez  [1]  derived  the  equivalent  sources  and  traveling  wave  current  on 
the  wire  excited  by  a  small  aperture  using  both  mode-matching  and  reci¬ 
procity,  but  he  'did  not  consider  the  interaction  of  the  wire  on  the 
electromagnetic  field  in  the  aperture.  For  a  small  circular  aperture 
whose  diameter  is  much  smaller  than  the  shortest  distance  from  the  wire, 
Lee  and  Yang  [2]  obtained  the  same  sources  as  those  derived  by  Kajfez 
plus  a  lumped  element  equivalent  circuit.  The  traveling  current  on  the 
wire  can  be  obtained  from  their  equivalent  circuit.  When  the  restriction 
on  the  size  of  the  aperture  relative  to  the  distance  from  the  wire  is 
removed,  Lee  and  Yang  modified  the  sources  but  not  the  lumped  elements. 
However,  the  lumped  circuit  elements  depend  upon  the  interaction  between 
the  wire  and  the  aperture,  which  should  be  taken  into  account  when  the 
wire  gets  closer  to  the  aperture. 

In  this  report,  we  obtain  a  pair  of  integral  equations  for  the 
problem  of  a  conducting  wire  passing  by  an  aperture  using  the  equivalence 
principle  [3].  These  equations  are  reduced  to  matrix  form  by  means  of 
the  moment  method  [4].  We  next  specialize  the  equations  to  the  problem 
of  an  arbitrarily-shaped  small  aperture  and  a  thin  inf initely-long  wire. 
The  electric  current  expansion  functions  on  the  wire  are  two  outward 
traveling  currents  (TEM)  on  each  semi-infinite  half  of  the  wire,  plus 
pulse  functions  representing  evanescent  currents  (higher-order  modes) 
in  a  finite  region  near  the  aperture.  The  magnetic  current  expansion 


j 
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functions  on  the  aperture  are  dipoles  located  at  the  center  of  the  aperture. 
For  simplicity,  a  plane  wave  excitation  from  the  region  opposite  the 
wire  is  assumed.  Finally,  an  equivalent  circuit  for  the  TEM  mode,  includ¬ 
ing  sources  and  lumped  elements,  is  derived.  We  obtain  both  the  outward 
traveling  currents  plus  evanescent  current  on  the  wire,  and  the  equivalent 
dipole  current  elements  of  the  small  aperture. 

II.  PROBLEM  FORMULATION 

The  general  problem  configuration  is  given  in  Fig.  2.1,  which  shows 
an  arbitrary  aperture  (labeled  A)  in  an  infinite  conducting  screen  separat¬ 
ing  regions  a  and  b.  In  region  a  (y < 0) ,  the  impressed  sources  J^a  and  M*a 
produce  short-circuited  fieLds  (fields  radiating  in  the  presence  of  the  com¬ 
plete  screen)  E13  and  H13.  In  region  b  (y  0)  ,  the  impressed  sources  J**5 
and  produce  short-circuited  fields  aruj  H^.  There  is  also  a  con¬ 
ducting  wire  labeled  B  in  region  b.  The  loss-free  homogeneous  material  is 
characterized  by  ;i  ,  in  region  a,  and  by  ;  ,  r  in  region  b.  The  problem 

3  cl  n  D 

is  to  find  the  currents  on  the  wire  and  in  the  aperture,  and  the  parameters 
of  the  equivalent  network  seen  by  the  TF.M  traveling  wave  on  the  wire. 

Following  the  development  of  [51,  we  treat  the  problem  as  a 
boundary  value  problem  for  which  tangential  components  of  both  J|  and  H 
are  continuous  across  the  aperture  and  the  tangential  component  of  E  on 
the  surface  of  the  wire  vanishes,  and  then  use  the  equivalence  principle 
[3].  In  region  b,  the  field  is  expressed  in  terms  of  the  equivalent 
magnetic  current  M  =  n  *  E  over  the  aperture  and  the  current  J_  on  the  wire. 
Here  M  corresponds  to  -M  in  [5],  n  is  the  unit  vector  normal  to  A  and 
pointing  to  B,  and  J5  is  the  electric  field  in  A  of  the  original  problem. 
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Fig.  2.1.  The  general  problem. 
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The  field  in  region  a  is  that  due  to  -M  plus  the  impressed  field.  To 

compute  M  and  J,  we  use  method  of  moments  [4]  and  expand  M  and  .J  as 
NA  MB 

M  =  J  V  M  and  J  =  y  I  J  ,  where  NA  and  NB  are  the  total  number  of 

—  u  n  — n  —  u  n  — n 


n=l 


n=l 


expansion  functions  over  A  and  B,  respectively,  and  V  ,  I  are  the  coef- 

n  n 

ficients  to  be  evaluated.  A  pair  of  matrix  equations  are  then  obtained, 


where 


fva 

V)  ►  ■>  ->  x 

+  Y  ]V  +  [T]  I  =  I 

(2-1) 

[T] V  +  fZ]I  =  V1 

(2-2) 

a 

a 

FYbl 

[■'"  NA  >■  NA 

(2-3) 

[T] 

-t(^n)  A^NA  x  NB 

(2-4) 

rf  i 

^t  B^NB x NA 

(2-5) 

[Z] 

f<~  ^m’  -t(^n)  B^NBXNB 

(2-6) 

I1 

=  ['M  ,  Hlb  -  H13  J...  v  . 

1  — m  -t  -t  A  NA  -  1 

(2-7) 

v1 

=  ^'4’  -t  B^NB  x  1 

(2-8) 

-V 

V 

^Vna  •  1 

(2-9) 

-> 

I 

f  IiJnB  >  1 

(2-10) 

Here,  M  and  J  are  the  testing  functions  for  M  and  J,  respectively. 

— m  - 1 m 

Subscripts  A  and  B  denote  symmetric  products  over  A  and  B,  respectively. 
Subscript  t  denotes  the  component  tangential  to  the  surface  (A  or  B) 
where  the  product  is  applied.  Superscript  a  denotes  region  a,  and 
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superscript  b  denotes  region  b.  (H(M  ),  E (M  )}  is  the  field  due  to 

- n - n 

M  ,  and  (H(J  ),  E(J  )}  is  the  field  due  to  J  •  Note  that  all  fields 
— n  - n - n  — n 

are  computed  with  the  aperture  shorted. 

If  a  Galerkin  solution  is  used  (i.e.  M  =  M  ,  J  =  J  ),  and 

— m  — m  — m  — m 

regions  a  and  b  of  Fig.  2-1  are  filled  with  the  same  material,  we  can 
use  reciprocity  [3,  sec.  3-8]  and  image  theory  [3,  sec.  3-4]  to  reduce 
(2-3)  -  (2-8)  to 


,b 

.a 

Y 

Y  =  2- -M 

mn 

nm  — m 

<-M  , 

H  (J  ,  .!')■• 

— m 

— t  — n  — n  i 

T  =  2  < -J  ,  E  (M  )>  =  -  T 

mn  — m  — t  — n  B  nm 


(2-11) 

(2-12) 

(2-13) 


Z 

z 

•  -J  , 

,  E  (J  , 

•J ' )  '• 

(2- 

-14) 

mn 

nm 

— i m 

— t  — n 

— n  B 

I1  = 

2->l  , 

Hiob 

-  Hioa • 

A 

(2- 

-15) 

m 

— i tn 

-t 

— t 

V1  = 

2'  J  , 

Ei°b. 

B 

(2- 

-16) 

m 

— t 

image 

of  J 
— n 

and 

located 

in  region  a. 

All  the  fields 

are  evaluated  with  the  conducting  plane  removed. 


III.  PROBLEM  SPECIALIZATION 

In  this  section,  we  specialize  the  problem  to  a  thin  infinitely- 
long  conducting  wire  of  radius  r  passing  by  a  small  arbitrarily-shaped 


aperture  centered  at  the  coordinate  origin.  Space  is  filled  with  loss- 
free  homogeneous  material  of  permeability  u  and  permittivity  e  on  both 
sides  of  the  screen.  The  impressed  field  in  region  a  is  assumed  to  be 
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a  plane  wave.  For  simplicity,  we  consider  horizontal  polarization  for 
the  incident  plane  wave.  The  incident  fields  with  the  screen  removed 
are 


Hi°a 


(3-1) 


ioa  -ikz 

F.  =  u  e 

-y 


(3-2) 


where  u  and  u  are  the  unit  vectors  in  the  x  and  the  y  directions, 

—x  — y 

respectively.  ”  =  denotes  the  intrinsic  impedance  of  free  space, 

and  k  =  2tt / >  where  '  is  the  wavelength.  There  is  no  incident  wave  from 
region  b.  The  configuration  of  the  specified  problem  is  shown  in  Fig.  3.1. 

A  Galerkin  solution  is  used  to  solve  the  problem,  and  we  specify  the 
magnetic  current  expansion  functions  in  the  aperture  and  electric  current 
expansion  functions  on  the  wire  as  follows. 


(A)  Magnetic  Current  Expansions  in  the  Small  Aperture 

As  described  in  If],  the  field  radiated  by  a  small  aperture  is  the 
combination  of  that  from  a  magnetic  dipole  tangential  to  the  screen, 
and  an  electric  dipole  p^  normal  to  the  screen,  both  located  at  the  center 
of  the  aperture  [7],  [8],  [9].  The  magnetic  current  Mon  the  aperture  can 
then  be  expressed  by  a  three-term  expansion, 

3 

M  =  )  V  M  (3-3) 

—  ,  n  — n 

n=l 

where  and  are  quasi-static  magnetic  currents  which  produces  the 
effects  of  unit  magnetic  dipoles  Kp  =  t^  and  =  t^,  respectively. 

Here  and  t ^  are  unit  vectors  tangential  to  the  screen  chosen  to 


I 

I 
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t 


Fig.  3.1.  A  thin  inf initely-long  wire  passing 
by  a  small  aperture. 
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diagonalize  the  magnetic  polarizability  *  .  For  simplicity  they  are 
assumed  in  the  x  and  the  z  directions.  is  a  quasi-static  electric 

current  producing  the  effect  of  an  electric  dipole  17  =  -  jur  n  .  Here 
n  is  a  unit  vector  normal  to  the  screen  and  is  chosen  in  the  y  direction. 
These  dipoles  are  located  at  center  0  of  the  aperture. 


(B)  Electric  Current  Expansions  on  the  Wire 


We  make  the  following  approximations  for  the  current  on  a  thin 
wire:  (a)  The  current  flows  only  in  the  axial  direction  of  the  wire, 
and  (b)  depends  only  on  the  axial  length  variable.  We  next  note  that 
the  wire  and  the  complete  screen  form  a  transmission  line.  We  there¬ 
fore  expect  that  the  aperture  excites  two  outward  traveling  (TEM)  waves 
propagating  in  the  -z  and  the  +z  directions  without  attenuation,  plus 
evanescent  waves  existing  near  the  aperture.  It  is  convenient  to  divide 
a  section  of  the  wire  near  the  aperture  into  NB-2  equally-spaced  sub¬ 
sections  of  length  2\  .  Pulse  functions  over  each  subsection  are  used. 

B 

The  expansion  functions  for  the  current  J  on  the  wire  are  then 


it " 

u 

—7. 

e.lkz 

z  :  0 

(3-4) 

— NB_ 

U 

— Z 

-  jkz 

e  •'  , 

0  •- 

z  <  * 

(3-5) 

J  = 
— n 

u  ^ 

P(z-z  )  , 

n 

n  = 

2,3,... 

,  NB-1 

(3-6) 

P(z-z  ) 
n 


1  for  z  _<  z  _<  z+ 
n  n 


(3-7) 


0  elsewhere 
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where  WL  is  the  length  of  the  section  over  which  the  evanescent  current 

is  assumed  to  exist,  z  ,  z  and  z,  are  the  z-eoordinates  of  the  start- 

—  n  + 

n  n  + 

ing  point  n,  midpoint  n  and  termination  n  of  the  nth  subsection,  and  the 

,  .  wl. 

numbering  starts  at  z  =  -  . 


IV.  EVALUATION  OF  MATRICES 

In  this  section  we  evaluate  t  ho  matrices  defined  in  (2-11)  -  (2-16) 
for  the  problem  of  Fig.  1.1. 


(A)  Impedance  Matrix  [Z] 

The  tangential  component  of  electric  field  on  the  surface  of  the 

wire  due  to  current  expansion  .1  at  the  rectangular  coordinate  (0,  d,  z’) 

— n 

and  its  image  I '  =  -  .1  at  (0,  -d,  z'l  are  expressed  in  terms  of  potential 
— n  —a 

integrals  as  (4,  f-.'qs.  (4-6)  -  (4-9)] 


-jkr  -jkr' 


F.  (J  -  j—-  d  (r’)  (- - 

z  — n  — n  4  n  r 


rr-)dz' 


'  d.l  (•/.')  -Jkr  -ikr ’ 
■  1  n  ,o  e 


”.)  •  ' z  j 

B 


1 — )  dz’ 


(4-1) 


,  n  =  1,2,3, - NB 


Here 


=  JoT 


•>  r> 

+  rg 


(4-2) 


is  the  distance  between  a  field  point  on  the  surface  of  the  wire  and  a 

source  point  of  .1  (z')  on  the  axis  of  the  wire,  and 
n 


z')Z 


+  (2d)‘ 


(4-3) 
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is  the  distance  between  the  field  point  and  the  image  of  the  source 
point . 

We  now  evaluate  the  elements  of  [Z]  matrix  defined  in  (2-14).  For 

m  4  1,  NB,  and  n  4  1,  NB,  it  is  the  [Z]  matrix  of  a  wire  of  length 

wL  +  2  a  above  a  ground  plane.  We  use  pulse  functions  for  both  J  (z') 
B  n 

3J  (z’> 

and  — ?— —  in  (4-1)  and  point  matching  for  testing  in  (2-14),  and  ex- 
tend  the  derivations  of  [4,  chap.  4]  to  include  the  due  to  the  image 
of  the  wire.  The  results  are 


Z  =  j.'J  V  '■  [  f  (  z  ,  z  ,  r  )  -  v(z  ,  z  ,  2d)] 
mn  mn  m  n  B  mn 


+  —  f.(z+,  7.+  ,  rB) 
J  m  n 


V(z+»  z+,  2d) 
m  n 


-  i(z+,  z_,  rfi)  +  w(z+,  z_,  2d) 
m  n  m  n 


-  ;(z_,  z+,  r  )  +  y(z_,  z+,  2d) 


m  n 


+  l(z_,  z_,  rR)  -  Mz_*  z_,  2d)] 
m  n  m  n 


where 


•i  (z  ,  z  ,  ‘ )  =  . 
m  n  4  i, 


1  f 


-)k 


2  2 

(z  -z’r  + 

•m  m 


dz’ 


n  A'  J(z  -  z’)2  +  p2 

n  ^  m 


(4-4) 


(4-5) 


has  been  evaluated  in  f 10 ] . 

For  m  =  1,  2,...,  NB  and  n  =  1,  NB,  from  (2-14)  we  have 


Z  =  -  J  ,  F.  (T  ,  J’)  ' 
m,  1  m  z  — I  —I  B 


(4-6) 


Zm,NB  ~  Jm’  V-NB’  ^NB)  B 


(4-7) 
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By  (11,  Eqs.  (A-9),  (A-12)]  and  (4-1),  we  have 


E 

z 


r  /JJSa 

llr  J  ‘  2 

o  r 

o 


r 

o 


~jkro 


-  (^k  - 
r' 


jkz 

9 


.3 


)  e 


-Ikr' 

o 


j  f  ,  jk  jkz  z 
__  [(--  +  —  +  -3)  e 


-ikr 


_  ,ih  +  1^4  + 


,3 


-jkr' 

)  e  °] 


where 


r 


o 


+  rl 


J->  o 

7.-  +  (2d)" 


However,  we  note  that  the  charge  densities  of  and  are 


dz 


=  ,1k  e 


ikz 


(z) 


z  <  0 


(4-8) 


(4-9) 

(4-10) 

(4-11) 

(4-12) 


_  .  . 

*  -  jk  e“,kZ  +  :(z)  0  <  z  <  ”  (4-13) 

dz 

and  E  due  to  these  point  charges  -  •  (0  )  and  (0+)  are  exactly  the 
— z 

-2  -3  -2  -3 

terms  involving  r  ,  r  ,  rf  and  r'  in  (4-8)  and  (4-9)  .  In  order 

o  o  o  o 

to  assure  the  continuity  of  the  total  current  (outward  traveling  current 

plus  evanescent  current)  on  the  wire,  we  first  replace  <?(z)  in  (4-12)  and 

P  (z) 

(4-13)  by  a  pulse  function  over  a  small  region  centered  at  z  =  0. 


d 
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Here  A.C^  is  chosen,  for  simplicity,  as  the  same  length  as  each  evanescent- 
current-subsection.  We  next  replace  the  E  due  to  charge  density  6(z)  by 

p  (z\ 

the  JS  due  to  charge  density  - - .  Thus  (4-8)  and  (4-9)  become 

0  -jkr  -jkr ' 


±i>  -  4^ 


-k  e. 


-) 


1 


e-Jkr  "jkr' 

AL 


0 

-■jkr  -jkr’ 

V^B’  ^  =  4^  <V^  - 


(4-14) 


4  i : 


1  ■  r  e"jkr  e"jkr’ 

'  v;o  z  r 


T— )  dz’  (4-15) 


where  r,  r',  r  and  r'  are  defined  in  (4-2),  (4-1),  (4-10)  and  (4-11), 
o  o 

respectively. 

Substituting  (4-14)  into  (4-6),  and  using  the  integration  formula 
[11,  Eq.  (37)],  we  obtain 


Z,,  =  —  ['.(0,  0,  r  )  -  !.T 0 ,  0,  2d)  1  (4-16) 

l  l  j  •  D 

*  FT  W>  ♦  ei<2ilt>  +  1 

[S.(rBk)  -  SjOdWP  --i- 

fa (o,  0,  r  )  -  ;(0,  0,  2d)]  (4-17) 

D 

k  A  • 

Z  .  =  - ^  [v(7.  ,  0,  r  )  -  j(z  ,  0,  2d)] 

ml  m  B  m 


(continued  on  next  pa£e) 
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i 


I 


+  W(z+.  0,  rB)  -  ^(z+,  0,  2d) 

J  m  m 

-  'M z_,  0,  rfi)  +  vjj<z_ ,  0,  2d)]  (4-18) 

m  m 

for  m  =  2,  3,  ....  NB-1  . 

Here  \p  is  defined  in  (4-5),  and  and  are  cosine  and  sine  integrals, 
respectively. 

Substituting  (4—1 5)  into  (4-7),  and  using  [11,  Eq.  (37)],  we  obtain 


ZNB,NB 

ZM,NB 


k  A.' 

m 

ojr 


[l  (z  ,  0, 


m 


V 


-  JJZ  [Hz 

-  h(z  ,  0 , 

m 

for  m  =  2,  3 


V 


rB)+0(z- 

m 

, . . . ,  NB-1 


-  ty(z  .  0,  2d)  ] 
m 

<Hz  ,  0,  2d) 

T 

m 

o,  2d)] 


(4-19) 


(4-20) 


Since  a  Galerkin  solution  is  used,  the  rest  of  the  impedance  elements  can 

be  obtained  using  Z  =  Z  for  m,  n  =  1,2,...,NB. 
mn  nm 

(B)  Source  Vectors  V*,  I*  and  Admittance  Matrix  [Y3  +  Y _ ] 

Since  =  £  in  (2-16),  it  is  obvious  that 

V1  =  0  (4-21) 

As  developed  in  [6],  the  source  vector  l1  and  admittance  matrix 
[Y3  +  yk]  for  a  small  aperture  in  Fig.  2.1  can  be  obtained  as  follows. 

We  first  interpret  as  an  electric  dipole  lo_  =  ).•  ^n  if  it  is  in 


I 
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I 


region  a  and  It  =  -ju<_  n  if  it  is  in  region  b,  while  M,  and  M«  are  the 

d—  —I  —  l 

same  as  those  defined  in  Section  III.  Applications  of  reciprocity 
[3,  Sec.  3-8]  to  (2-7)  and  the  use  of  the  definitions  of  M^,  and 
give  the  source  vector  T 

T1  =  H^fO)  -  (0)  (4-22) 


T2  =  H2b(^  '  H2S(-) 


(4-23) 


(rbE3b»>  -  raEia<WI 


(4-24) 


Here,  and  denote  the  t^  and  the  components  of  H,  and  denotes 

the  n  component  of  E.  All  the  fields  are  evaluated  at  the  center  ()  of  the 
aperture . 

The  imaginary  part  of  fY3  +  Y^]  can  be  evaluated  by  the  extension  of  the 

Bethe-hole  theory  [7]  to  allow  for  contrasting  media  [12].  Thus,  the  field 

due  to  -M  in  region  a  of  Fig.  2.3  can  be  approximated  by  the  combination 

of  that  of  a  magnetic  dipole  with  current  element  -K  plus  that  of  an 

electric  dipole  with  current  element  -I ■ .  The  field  due  to  M  in  region  b 

of  Fig.  2.1  can  be  approximated  by  the  combination  of  that  of  a  magnetic 

dipole  with  current  element  K  '  =  plus  that  of  an  electric  dipole  with 

current  element  T  •  '  In  addition,  the  relationships  among  dipole 

a 

moments,  incident  fields  and  these  current  elements  are 

=  ’-a  B™  (4-25) 


l  P 


(4-26) 
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2u. 

- T -  {-  t,  1 

U  +  U,  —1  ml 
a  b 


[Hja(0)  -  Hjb(0)] 


-  t2  -  H2b(0)]  • 

(4-27) 

2-  •  •  v 

P  =  - ~~  n  a  [  E^a(0)  -  c  E^b(0)l 

— e  +  C  —  e  a  3  —  b  3  — 

(4-28) 

a  b 


Here  a  ,  and  a  „  are  the  t,  and  t_  components  of  magnetic  polariza- 
ml  m2  —1  —2 

bility,  and  a  the  electric  polarizability.  The  relationships  between  the 
coefficient  vector  V  of  M  and  these  current  elements  K?,  I  >'  can  be  obtained 
from  the  definitions  of  ,  M„  and  .  They  are 

Kl  =  t_  +  V2  t2  (4-29) 

_K  =  -.1  ■  a  V  n  (4-30) 


->■  ►  b 

By  (2-1)  with  T  =  0  ([Y  +  Y  ]  can  be  obtained  without  the  presence  of 

the  wire)  and  (4-22)  -  (4-30),  we  have 

u  + 


[va  +  ybi 


2J-; 


‘a “b  ml 


0 


+  ;i 

a  b 

2j.;,  \i  t  ~ 
a  b  m2 


j*'  C  +  o 

a  b 

2a 


(4-31) 


The  admittance  matrix  in  (4-31)  came  out  purely  imaginary  because  no 
radiation  terms  were  included  in  the  equivalent  dipoles  of  the  small 
aperture  in  (4-27)  and  (4-28),  that  is,  only  the  first  term  of  a  fre¬ 
quency  expansion  for  [Ya  +  Yl  j  is  evaluated.  When  there  is  an  object 
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V 


I 


« 


near  the  aperture,  as  in  our  original  problem,  the  first  term  of  expansion 
3  b 

for  [Y  +  Y  )  can  be  cancelled  by  the  interaction  between  the  object  and  the 

aperture.  Then  the  second  term  of  frequency  expansion  for  [Ya  +  Yb] ,  which 

is  its  real  part,  becomes  important. 

3  b 

The  real  part  of  fY  +  Y  ]  can  be  obtained  from  the  power  radiated  by  M. 

Note  that  M  are  real  magnetic  currents  (M  equivalent  to  a  magnetic  current 
— T)  —  J 

loop  KS  =  — —  [3,  p.  1351),  and  for  M  real,  we  have  Re(Ya  +  Yb  )  equal  to 
.t.  —n  nn  nn 

b 

the  power  radiated  bv  M  in  regions  a  and  b.  Using  the  formula  for  the  power 

— n 

radiated  by  an  electric  dipole  [3,  Kq .  (2-116)1,  duality  [3,  sec.  3-2]  and 

definitions  of  >1  .  we  obtain 
— n 


2  2 
k  k 

Re  ( Y  .j1 1 )  +  Re  (Yb  ]  )  =  -77  (~  +  7—) 


11 


(4-32) 


2  2 

Re(Ya2)  +  Re(Y22)  =  p  (7^  +  7^) 

a  b 


4  4 

Re(Ya  )  +  Re(Yb  )  =  jr  (~  +  7^) 

a  b 


(4-33) 


(4-34) 


2 ,,  2  • 

where  ^  =  1/ ^  =  —  ,  kfe  =  —  and  ^  . 

“a  a  b  » 

(4-34)  we  obtain  the  admittance  matrix  for  a  small  aperture  in  Fig.  2.1. 

For  the  problem  of  Fig.  3.1,  we  have  ,  =  {,  ’  a  b 

Li  *  u  =  .  Thus ,  the  source  vector  I1  can  be  obtained  by  (4-22)-(4-24) 

a  b 

with  Hlb  =  0,  KLb  =  0  and  Hia,  F.La  defined  in  (3-1)  and  (3-2).  The  re- 


.  Combining  (4-31)  - 


=  c ,  \  =  X  =  A  and 


suit  is 


"il 

— 

2 

T, 

1 

?*> 

ij 

0 

.'i. 

:2.i  •• 

(4-35) 
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The  admittance  matrix  of  (4-31 )-(4-34)  can  be  reduced  to 


8 "  +  1 


3<,A2  jU"':  "ml 


(Y3  +  Yb] 


8  +  1 


3n 


\2  j(,,y 


m2 


12_  + 

-  ,  4  oi 
3  HA  e 


(4-36) 


(C)  Coupling  Matrices  [ T 1  and  [ T ] 

Since  a  Galerkin  solution  is  used,  [T]  can  be  obtained  in  terms 
of  [T]  according  to  (2-13).  Therefore,  we  only  have  to  evaluate  [T] . 
Application  of  reciprocity  to  the  right-hand  side  of  (2-12)  gives 


T 

mn 


ff 

.1 

jJ  ~T1 

B 


E(2M 

—  — T) 


)ds 


(4-37) 


Once  again,  application  of  reciprocity  to  the  right-hand  side  of  (4-37) 
and  use  of  the  impulsive  nature  of  M  give 


T 

mn 


-H  (.1  ,  J'  ;  0) 
m  — n  — n  — 


E  ( J 
3  — r 


J’  ;  0) 


m  =  1,2 


m  =  3 


for  n  =  1,2,...,  N'B 


(4-38) 


Here  H  (J  ,  J':  0)  and  E„(J  ,  J';  0)  are  the  fields  due  to  J  and  its 
— m  — n  — n  —  3^n~n—  — n 

image  J'  evaluated  at  0. 

— n  - 

The  x  and  the  z  components  of  H  fields  due  to  Jj  and  at  the 

center  0  of  the  aperture  can  be  obtained  by  [11,  Eqs.  (A-6)  and  ( A— 11)] 


with  the  coordinate  of  ^  at  p  =  d,  <0  =  —  and  z  =  0  in  the  geometry 


system  of  [11,  Fig.  A— 1 ] . 


H  (J  ;  0)  =  H  (.T  ;  0) 

x  —  1  —  x  — NB  — 


=  -  sin  ^  H  (Jj  ;  0) 

—  1  kd 
e 

4"d 


■  °>  ■  V4b  >  ® 


The  v  component  of  JE  fields  due  to  can  be  obtained  by  using  [11, 


Eq.  (A-ll ) ] ,  giving 


E  (.T.  )  =  sin  C  E  (J.  ) 
v  —1  —1 


sinf 

.i  ‘ 


sin 


+  jM  _  -1kz 


2  "ikrn 

2  i  O 

-  -3]  c 


where  r 

o 


Again,  at  0  wc  have 


d,  0 


32 

2 


and  z  = 


(4-41)  becomes 


0) 


-  jkd 
•e 


d 


2 


Similarly,  using  [11,  Eq .  (A— 6 ) J ,  we  have 


Ev<4.B  ;  2>  ■ 


t  (Jj  ;  0) 


The  ^-component  of  H  due  to  J  for  n  =  2 , 3 , . . . ,NB-1  are 
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H  ,(J  )  = 
(p — n 


A  (J  ) 

z  — n 


1  •  f 


-jk  'l(z-z')2  +  n2 


J  2 

yj(z-z’)  +  ' 


dz* 


V>  ,  -  jkr 

— — —  (J_  +  iiS)  e  n 

4"  v3  2  ^ 

r  r“~ 

n  n 


(4-44) 


where  r 


=  a  (z-z  ) 2  + 


n  ' '  n 
evaluated  at  0  are 


Thus,  the  x  and  the  z  components  of  H(J  ) 

—  — n 


H  (.T  ;  0)  =  -  sin  1  H.  (  T,  ;  0) 

x  — n  —  .  — 1  — 


d  ,  -jkr 

n  ,  1  ]k  on 

“  (~r  c 

r 


on  on 


(4-45) 


H  (J  ;  0)  =  0 
z  — n  — 


(4-46) 


where 


r  = 

•  m  v  n 


+  d‘ 


(4-47) 


The  y  component  of  K  due  to  at  0  can  he  obtained  by  ui  mg 
sin  * 

E  =  — ; - 1  —  H,  evaluated  at  0,  giving 

V  J  •  '7.  . 


K  (d  ;  0)  =  7-rr—  {(■—  +  e 
v—n  -  a  .1  •  3  i 

+  + 

n  n 


•jkr 


-jkr 


<W 

r  r 


(4-48) 


for  n  =  2,3,..., 


NB-1  . 
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where 


We  next  note  that 


r+ 

n  n 


r  = 

n 


n 

~T 

aJz_ 


+  d 


+  d 


(4-49) 


(4-50) 


H  (.1  ,  .1'  ;  0)  =  2  H  (J  ;  0) 

x  — n  — n  —  x  — n  - 


(4-51) 


K  (J  ,  J'  ;  0)  =  2  E  (J  ;  0) 
v  — n  — n  —  v  — n  — 


(4-52) 


for  n  =  1  ,  2  , .  .  . ,  NB  . 

Therefore,  by  (4-381  -  (4-40),  (4-42),  (4-43),  (4-45).  (4-46),  (4-48), 
(4-51)  and  (4-52),  we  have 

-jkd 


T  =  T 
11  l.NB 


-e 


T  -  r  = 

31  3.NB 


-d 

-jkd 

e 

-> 

2-  .r 


•  u  . 

_n  /  J_  +  jJL 


In  2  •  '3 

r  r 
on  on 


)e 


-j  k  r 


on 


-jkr+ 

d  .  ,  1  jk.  n 

r3n  =  F  1  3  +  c 

r+  r+ 

n  n 

-jkr_ 

,  1  ,  jk.  n, 

-  (-3  +  F)  c  ]  , 

r  F 


(4-53) 

(4-54) 

(4-55) 


(4-56) 


for  n  =  2,3,...,  NB-1 


T-  =  0  ,  for  n=l  ,2 . NB 

2n 


(4-57) 


Finally,  we  note  that  for  the  chosen  incident  field,  substituion  of 
ij  =  0  in  (4-35)  and  =  0  in  (4-571  into  (2-1)  gives  =  0.  This  means 
that  there  is  no  equivalent  dipole  in  the  J. direction  for  the  incident 
field  specified  in  Fig.  3.1.  In  order  to  save  computation,  we  expand  M  as 
a  two-term  expansion  with  representing  Kj_  =  t_  and  H  representing 
1  <_  =  -  jwt  n,  and  remove  the  corresponding  columns  or  rows  due  to  Ki  =  _t 
in  the  matrices  I*,  [Y'1  +  Y^'l,  [T],  f T ]  and  V. 


EQUIVALENT  ClRCl’TT 


In  this  section,  an  equivalent  circuit  for  the  transmission  line  mode 
on  the  wire  at  ?.  =  0  is  derived.  Basically,  the  approach  is  similar  to  that 
by  Yang  and  Harrington  [111  who  dealt  with  the  problem  of  a  narrow  slot 
backed  by  an  inf  in i tel v-1 ong  wire.  Actually,  the  network  they  derived  is 


good  onlv  for  an  aperture-wire  problem  for  which  I.  =  I  .  In  the  follow- 
ing  we  generalize  the  network  such  that  it  can  he  applied  to  any  aperture- 
wire  problem.  It  is  known  that  an  equivalent  circuit  of  a  wire  passing 
by  an  arbitrary  aperture  is  a  two  port  network  as  shown  in  Fig.  5.1.  The 


parameters  [7/]  depend  onlv  on  the  geometry  of  the  problem.  They  can  be 


obtained  bv  removing  the  incident  fields  from  the  system  and  applying 
mathematically  arbitrary  excitations  to  the  transmission  line  formed  by 


e 

the  wire  and  the  aperture-perforated  conducting  plane.  The  sources 
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(A)  Network  Parameters  [ Z*-  ] 

To  obtain  the  elements  of  [  7.^  ]  ,  we  use  the  definitions 

(5-1) 

(5-2) 

where  i^ ,  i0,  and  v,  are  port  currents  and  voltages  at  ports  1  (z=0  ) 

and  2  (z=0+),  respectively,  with  reference  directions  shown  in  Fig.  5.2. 

To  represent  an  excitation,  we  mathematically  apply  a  TEM  current 
i  k  z 

wave  I  e  J  u  to  the  transmission  line.  The  corresponding  incident 
— z 

fields  from  region  b  of  Fig.  2.1  are  then 


vi  =  zn  \  4  zi 2  *2 


V2  Z21  ' 1  +  Z22  12 


Hib  *  T+  h  e"ikz 


^ib  +  -ikz 
E  =  V  e  e 


(5-3) 

(5-4) 


where  T+  and  V+  are  the  mode  current  and  voltage,  and  e  and  h  a  e 

transverse  electric  and  magnetic  fields  of  Tt  1  •  -ode .  *e  now  have  an 

+  - i kz  ,  ,  .  .  p  j  kz 

incident  wave  I  e  u  ,  and  two  outward  traveling  currents  I.  e  u 

— 7.  1  — z 

and  I  e  ^kz  u  excited  bv  the  aperture,  as  shovTi  in  Fig.  5.3.  I.  and 
NB  — z  i 

I  can  be  obtained  by  solving  (2-1)  and  (2-2)  with  the  incident  fields 
NB 

defined  in  (5-3)  and  (5-4). 

Combining  Fig.  5.2  and  Fig.  5.3,  we  obtain  the  port  currents  i^  and 
i2,  and  port  voltages  v^  and  v.,  at  ports  1  and  2,  respectively.  They  are 


1 


T+  +  1. 


a  +  V 


-  Z  (I  -  I.) 

o  1 


z  0  +  iMR) 

o  NB 


(5-5) 

(5-6) 

(5-7) 


V 


(5-8) 


Fig.  5.3.  A  TEM  incident  current  wave  I  e  u  produces  two 

ikz..  T  z-jk2.. 


outgoing  currents  1  e 


\  and  IN'B  e 
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where  Z  =  60  log  ( — )  is  the  characteristic  impedance  of  the  trans- 
°  rB 

0 

mission  line.  Since  fZ  ]  is  independent  of  the  excitation,  (5-1) 
and  (5-2)  hold  for  the  port  currents  and  voltages  in  (5-5)  -  (5-8). 
That  is 


=  Z 


11  ‘1 


+  Z 


12 


(5-9) 


=  Z 
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'  .  -,c 

X1  +  Z22 


(5-10) 


0 

To  solve  for  [Z  ],  we  need  two  more  equations  relating  the 

0 

port  currents  and  voltages  to  fZ  ]  in  addition  to  (5-9)  and  (5-10). 

—  i  k  z 

One  can  apply  a  TEM  current  wave  T  c>'  «£  to  the  wire.  Again, 

o  i  k  z 

the  aperture  excites  two  outgoing  travel inq  currents  I  eJ  u 

1  z 

—  i  k  2 

and  I  e  J  u  ,  as  shown  in  Fie.  5.1*.  The  port  currents  and 

—7. 

voltages  for  this  case  can  he  obtained  hv  combining  Fig.  5.2  and 
Fig.  5.4,  giving 


41 = 

■' + 

(5-11) 

*  2  ~ 

-  <r  +  V 

(5-12) 

vi  = 

- V" +  h> 

(5-13) 

V2  = 

z  (L,  -  r> 

i'  NB 

(5-14) 
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Similarly  to  (5-9)  and  (5-10),  we  have  a  pair  of  equations  re- 

e 

lating  these  port  currents  and  voltages  to  [Z  ].  They  are 


~  „e  i  ,  _e  ? 

V1  =  Z11  X1  +  Z12  X2 


(5-15) 


'  e  3  e  * 

V2  =  Z21  *1  +  Z22  X2 


(5-16) 


Thus,  by  solving  (5-9),  (5-10),  (5-15),  and  (5-16)  with  port  currents  and 
voltages  defined  in  (5-5)  -  (5-8)  and  (5-11)  -  (5-14),  we  can  ob¬ 
tain  f  ze  ] . 

However,  we  note  that  if 


1  =  -  I 


(5-17) 


1 1  "  lNB 


(5-18) 


'xb  =  '  ! 


(5-19) 


This  can  be  verified  bv  introducing  a  new  coordinate  svstem  where  the  z 

—  •”  i  k  Z 

coordinate  is  replaced  bv  -?.  in  Fig.  5.4.  Thus,  we  have  I  e  (-u ) 

producing  1,  e  (-u  )  and  1  e'^z  (_u  ).  Comparing  these  currents 

in  the  new  coordinate  system  with  those  in  Fig.  5.3,  we  obtain  (5-17)  - 
(5-19).  Substituting  (5-17)  -  (5-19)  into  (5-11)  -  (5-14),  we  obtain 


*1  =  x2 


(5-20) 


X2  X\ 


(5-21) 


V1  “  V2 


(5-22) 


V2  =  V1 


(5-23) 


Substituting  (5-20)  -  (5-23)  into  (5-15)  and  (5-16),  we  obtain 
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V2  ZU  l2  +  Z12  Xl 


(5-24) 


V1  =  Z21  j2  +  Z22 


(5-25) 


By  (5-9),  (5-10),  (5-24)  and  (5-25)  we  have 


e  „e 
7.  =  7 

11  21 


Vlll  V212 


i]  -  x-2 


(5-26) 


/■  ,  y  .Vu'.Vi 
12  21  ;2  ;2 
31  '  l2 


(5-27) 


We  therefore  obtain  the  network  cleme-nts  hv  simply  applying  a  TEM  current 

wave  T  e  J  u  and  solving  for  the  outward  traveling  currents  or.  the 
—7. 


Furthermore,  (5-261  and  (5-27)  show  that  the  network  is  symmetric 
(Z®.|  =  Z^2)  and  reciprocal  (7.^  =  7^),  and  therefore  it  is  equivalent  to  a 
symmetric  Tee  network  as  shown  in  Fif».  5.5.  The  network  elements  are 


7°  =  7° 

1  11 


-  (I1  +  1NB)  Z0 

21+  +  h  +  ‘kb 


(5-28) 


Z6  =  7e 


*  'sb>  zo _ 

”,  -  V<:1++  h  +  V 


(5-29) 
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(B)  Sources  and  V® 


If  both  ports  of  Fig.  5.5  are  terminated  with  matching  loads,  as 

shown  in  Fig.  5.6,  the  problem  corresponds  to  that  of  an  inf initely-long 

wire  passing  by  an  aperture.  This  is  our  original  problem.  In  this  case, 

we  can  solve  matrix  equations  (2-1)  and  (2-2)  to  obtain  the  outward  travel- 

ing  currents  I,  eJ  '  u  and  1„.„  e  u  .  The  loop  currents  i,  and  i„  in 

1  — z  NB  — z  12 

Fig.  5.6  are  equal  to  T  and  -I  ,  respectively.  Thus,  Fig.  5.6  gives 

I  N  B 

the  sources 


V?  =  -  (Z0  +  7M  +  7-2)!l  +  Z?  lNB 


(5-30) 


■’2  =  (Z0  +  Z1  +  V 


e  e 

1)1  -  Z  I 

?’  NB  2  1 


(5-31) 


Note  that  Zj,  Z^,  and  in  (5-28)  -  (5-31)  are  proportional  to  Zfl. 


We  therefore  can  normalize  them  by  Z^.  This  completes  the  derivations 
of  the  equivalent  network  for  the  general  problem  in  Fig.  2.1. 

Because  of  the  1  inoar  it\  of  the  operators  in  (2-1)  and  (2-2),  1^ 

■s  4*  + 

and  T  _  are  proportional  to  I  .  Thus,  we  can  choose  I  =  1  and  reduce 
NB 

(5-28)  and  (5-29)  to 


c  "  °'l  +  iNB)  Zf 


Z1  =~ 


(5-32) 


2  +  !1  +  !nb 


2  (1  +  W  zo 


(I1  -  rNB)(2  +  !1  +  V 


(5-33) 


In  a  summary,  the  network  is  a  symmetric  Tee  network  (Fig.  5.5), 
and  the  elements  defined  in  (5-32)  and  (5-33)  can  be  obtained  by  removing 
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—  I  Ic  z 

the  incident  fields  from  the  system  and  applying  a  TKM  current  e 

to  the  wire.  Final lv,  wc  note  that  for  the  symmetric  case  Uj 

I  _  i  2°  in  (5-331  is  infinite  and  the  equivalent  network  in  Fig.  5.5 

1  NB  ’  2 

can  be  reduced  to  [11.  Fig.  8]  with  the  element  and  source  defined  in 
[11,  Eq.  (99)  and  (100)] . 

For  the  small  aperture  problem  ot  Fig.  3.1,  the  source  vectors 
fi  and  v1  corresponding  to  the  excitation  !+  e",k?'  mathematically  applied 

to  the  wire  are  derived  in  the  following,  wav. 

The  transverse  electric  and  magnetic  fields  of  TF.M  mode  at  center  0 

of  the  small  aperture  are  [1] 


h  (0)  '  - 


(5-34) 


e  ( 0  '  =  h(0)  •  u  ^ 


(  >-33) 


Bv  (4-2  21  -  (4-24)  with  "  '  "•  1 


* ,l  I1,1  =  0  and  (4-3),  (5-4),  (5-3-t) 


and  (5-35),  wv  have 


V  =  0 


(5-3b) 


(5-37) 


-i  1; 


In  (5-37),  V+  =  ZQI+  and  T+  =  1  are  used.  Note  that  12  ~  0  in  37) 

and  T  =  0  in  (4-57).  Hence,  as  stated  before,  we  can  reduce  N  to  a 
1  2n 

two  term  expansion  and  remove  the  corresponding  rows  or  columns  in 
[Y3  +  Yb],  IT],  IT],  T1  and  V. 


VI.  NUMERICAL  RESULTS  AND  DISCUSSION 

A  computer  propram  was  developed  to  calculate  the  current 
distribution  on  an  inf  ini  tel v-long  wire,  the  current  elements  of  the 
equivalent  dipoles  for  an  arbitrarily-shaped  small  aperture  whose 
polarizabilities  are  given,  and  an  equivalent  network,  for  the  trans¬ 
mission  line  mode.  Fortunatelv,  analytical  expressions  of  polarizabili¬ 
ties  for  circular  aperture  [71  and  elliptical  apertures  [9]  exist,  and 
numerical  solutions  for  any  other  shaped  apertures  have  been  developed 
[121,  [13],  [14].  In  this  section,  we  present  the  computational  results 
for  the  problem  of  a  small  circular  aperture  of  radius  r^  backed  by  a 

thin  infinitelv-lonc  wire  of  radius  rn  at  v  =  d. 

n 

In  Fir,.  6.1,  we  show  the  distribution  of  the  total  and  the 

evanescent  currents  on  the  wire,  for  the  case  of  r,  =  r_  =  0.001 >  and 

A  a 

d  =  0.1'.  The  real  and  the  imaginary  parts  of  total  current  are  shown  in 
Fig.  6.1(a)  and  (M,  respect  ive  1  v .  The  real  and  the  imaginarv  parts  of 
evanescent  current  are  shown  in  Fie.  6.1(c)  and  (d),  respectively.  It  is 
seen  that  the  evanescent  current  rapidly  drops  to  zero  beyond  a  small 
region  near  the  aperture.  Furthermore ,  both  the  outward  traveling  cur¬ 
rent  and  the  evanescent  current  ire  discontinuous  at  the  location  of  the 
aperture  (z  =  0),  but  the  total  current  is  continuous  everywhere  along 
the  wire.  This  shows  that  the  number  NB  of  the  current  expansions  on 
the  wire  must  be  even,  and  the  continuity  of  the  total  current  is  assured 
by  replacing  the  point  charge  ' (z)  bv  a  pulse  function  as  stated  before. 
In  fact,  the  discontinuities  of  the  evanescent  and  the  outgoing  traveling 
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I 


(Fig.  6.1.  The  current  distribution  on  the  wire  induced  by  a  small 

circular  aperture,  for  r^  =  r^  =  0.001'’  and  d  =  0.1\. 

(a)  Tlu  real  part  ot  total  current,  (b)  The  imaginary  part 

I  of  total  current.  <c>  The  real  part  of  evanescent  current. 

fd)  The  imapinarv  p  irt  oi  evanescent  current. 

I 

I 
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current  at  z  =  0  are  expected  from  the  results  in  Appendix  A,  since  =  I 

for  magnetic  dipole  problem  and  I.  =  -  I  for  electric  dipole  problem. 

In  Fig.  6.2,  we  illustrate  the  variations  of  the  equivalent  circuit 

parameters  when  the  distance  d  between  the  wire  and  the  aperture  is  changed, 

for  the  cases  r,  =  0.01  and  r  =  0.001  .  Figure  6.2(a)  shows  that  the 

imaginary  part  of  approaches  zero  when  *1  ■  r  .  The  real  part  is  much 

smaller  than  the  imaginary  part  and  can  be  neglected  when  d  ■  r  .  Figure 

—  A 

6.2(h)  shows  that  the  imaginary  part  of  Zf,  becomes  larger  when  d  is  larger, 
and  the  real  part  can  he  neglected  when  d  r^.  Figure  6.2(c)  shows  that 

C  t1 

Vj  is  the  complex  conjugate  of  y  (approx  imat  el  v ) . 

Before  we  compare  our  results  with  those  obtained  by  using  Kajfez's 

fl],  and  Lee  and  Vang's  [21  formulations,  wt  show  their  equivalent  networks 

in  Fig.  6.3.  In  their  formul at  ions,  d  is  the  distance  of  the  wire  to  the 

screen,  the  distance  of  the  ••■ire  to  the  center  0  of  th.e  aperture  and 

Rq  =  d  for  the  problem  considered  in  Li.;.  3.1.  If  both  terminals  are 

loaded  with  Z  ,  the  problem  is  equivalent  to  an  inf  in  i  to  1 v- long  wire  and 
o 

the  currents  anil  i ,,  at  z  =  0  md  u  in  lie.  h.3  correspond  to  I  ^  and 
-I  ,  respective!  v.  LV  therefore  can  -Main  and  i0  from  their  circuits 

and  compare  the  results  with  ou«-  I  .,nd  -I  .  Ihev  use  the  Bethe-hole 
theory  [71  to  approx ir at  e  the  equivalent  dipoles  of  the  small  aperture, 
but  we  have  additional  terms  to  account  for  the  radiation  neglected  in  the 
B«  the- hole  theorv.  It  is  i  n  t  ■  -r.  st  i  m  t  •  eo-.pare  the  current  elements  of 
equivalent  dipoles  using  both  •  -rmulas.  Kv  Bethe-hole  theorv,  the  current 


elements  of  equiva'ent  dinol* 


rc,-  i  -n  h  o'  hi-.  3.1  (v  0)  are 


nt  impedance  7.^  for  r  =  0.01  \  and 
(  =  1  meter,  E10'  =  1  volt/meter) 


I 


‘  30000 


20000 


9 


10000 


00 


Fifi.  6.2(h). 


The  equivalent  impedance  Y.\  for  r 

J  .  A 

r  =  0 . 00 1  •  .  (  •  =  1  meter.  F.U>a  = 


=  0.01 '  and 
1  volt/meter) 
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K;  =  -J  'mx  H^>  (6-1) 

1  ■  =  .)'«  !\!a  (0)  (6-2) 

The  corresponding  current  elements  of  our  solution  are 

K'  =  Vj  (6-3) 

1-  =  -j  V,  (6-4) 

Here  \'  and  V?  are  t  lie  coefficients  for  M  and  M  is  reduced  to  a  two-term 
expansion.  Comparing  (6-1)  and  (6-2)  with  (6-3)  and  (6-4),  we  can  observe 
that  the  equivalent  dipoles  hv  the  Rethe-hole  theory  depend  only  on  the 
incident  fields  and  tlu  dimensions  of  the  aperture,  but  the  equivalent 
dipoles  by  our  theory  depend  on  the  coup! ins  In-tween  the  wire  and  the 
aperture  ((!],  (Cl),  the  geometries  of  the  wire  (fZp  and  the  aperture 
([Y]),  and  incident  fields  ( 1 ' ) . 

We  next  compare  our  re.-u1'--  •■•ith  those  i-htained  hv  using  other  methods. 

In  Table  6.1,  we  compare  our  results  of  amplitudes  of  outward  traveling 

currents  on  the  wire  with  those  obtained  by  using  Kajfeg.'s,  and  Lee  and 

Yang's  formulas,  for  the  cases  r,  -  r  -  0.001  and  r,  =  0.01-,  r„  =  0.001'. 

A  ;>  A  d 

It  is  seen  that  our  results  agree  closelv  witli  the  cither  two  results  for 

the  cases  of  d  r  ^ ,  hut  disagree  with  the  other  two  results  f<ir  the  cases 

of  d  •  r  .  Since  all  three  methods  treat  a  small  aperture  as  dipoles, 

for  the  cases  of  d  r  such  approx imat ions  may  not  be  accurate.  There- 

A 

fore,  ,i  different  method  should  be-  used  for  these  cases. 

fable  6.2  shows  comparisons  between  our  results  of  current  elements 
(if  the  equivalent  dipoles  and  those  obtained  hv  using  the  Bethe-hole 


itudes 


N'o  correction  factor  is  used  in  Lee  and  Yang's  results. 
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Table  6.2.  The  current  elements  Ky_  and  _V_  of  equivalent  dipoles  in 
region  b  (y  0),  for  r  =  0.00]'.  N  is  a  constant  to  be 
divided  hv  K  or  f  . 


V1 

X 

d/1 
h ~ 

VI./  1 

0.001 

r-*- 

1 

O 

0.10 

J 

0.05 

1  .0 

62 

0.01 

0.5 
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o.oos 

O.l 

6J 
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io~4 

0.10 

7.0 

82 

0.05 

1  .0 

'02 

i 

0.02 

0.2 
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0.005 

0.05 

>  J 
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theorv,  for  the  cases  r  =  r  =  0.001 '  and  r  =  0.01',  r  =  0.001A.  It  is 

Ad  A  d 

seen  that  our  results  agree  (exactly)  with  Bethe's  results  for  the  cases 
when  d  r  ,  but  disagree  with  Bethe's  results  and  the  real  parts  of  the 

A 

current  elements  become  important  when  d  r  .  We  note  that  these  real  parts 

~  A 

correspond  to  the  radiation  from  the  equivalent  dipoles.  This  can  be  veri¬ 
fied  by  using  (2-1)  with  I  =  fi,  (4-15),  (4-16),  (0-3),  (6-4)  and  recalling 
that  the  real  part  of  [V''  +  V']  accounts  for  the  radiation  from  the  dipoles. 

From  the  comparisons  in  Tables  h.l  and  h.2,  we  can  conclude  that 
Kajfez's,  Lee  and  Yang's  and  our  solutions  are  good,  for  the  cases  of  d  •  r  . 

As  mentioned  prev  iousl v ,  all  the  three  solut  ions  may  not  be  accurate  for  the 
cases  of  d  r^.  However,  we  take  into  account  both  the  radiation  from  the 
equivalent  dipoles  and  the  interaction  of  the  wire  on  the  field  in  the  aper¬ 
ture.  Therefore,  wo  believe  111,0  our  solution  can  he  used  as  an  approxima¬ 
tion  for  the  cases  of  d  r  . 

For  the  problem  of  a  leaded  wire  pass  in  t  bv  a  small  aperture,  the  voltage 
and  current  of  TFM  wave  on  the  wire  can  be  obtained  by  applying,  typical  trans¬ 
mission  line  tormulat  ions  to  the  equivalent  circuit  for  transmission  line  mode. 
Further  work  on  solvino  tin.  tot  it  .  urtunt  on  the  loaded  wire  andtiie  current 
elements  of  equivalent  dipoles  of  the  small  aperture  can  be  done  using  a 
similar  approach  to  that  developed  in  [15].  A  general  solution  to  a  wire 
passing  by  art  arb  i  t rar i 1 v— si  zed  aperture  can  be  done  using  some  01  the  basic 
formulations  developed  here  and  aoplving  triangular  patching.  We  expect  L< 
report  the  results  in  the  m  nr  future. 
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APPENDIX  A 

DIPOLE  BACKED  BY  AN  INFINITE  WIRE 

We  here  consider  appl ications  of  (2-1)  and  (2-2)  when  the 
magnetic  current  in  the  aperture  is  a  given  magnetic  dipole  =  u 
or  an  electric  dipole  1=  j  u  ,  which  is  located  at  the  origin  of  the 

-  — v 

coordinate  and  hacked  by  a  z-directed  infinite  wire  of  radius  r^  at 
y  =  d.  For  this,  we  have  M  expanded  as  before  with  (V^  =  1,  =  0} 

for  the  magnetic  dipole  case  or  =  0,  V0  =  -  v}  for  the  electric 
dipole  case,  and  .1  on  the  wire  expanded  as  before.  Then  the  same 
analysis  as  that  used  for  a  small  aperture  can  be  applied.  Thus  (2-2) 

j 

with  V  =0  is  used  to  solve  for  I. 

Figure  A.l  shows  that  the  current  distribution  on  the  wire 
induced  by  the  magnetic  dipole  is  symmetric  about  z  =  0.  Figure  A. 2 
shows  that  the  current  distribut  ion  on  the  wire  induced  by  the  elec¬ 
tric  dipole  is  anti-svmmetric  about  z  =  0.  They  also  show  that  the 
evanescent  current  exists  within  a  small  region  on  the  wire  near  the 
dipole. 

Table  A.l  shows  that  our  results  agree  closely  with  those 
obtained  by  using  Kaj fez's  formulation  [11,  except  for  the  small 
imaginary  part  for  the  magnetic  dipole  case  and  the  small  real  part 
for  the  electric  dipole  case,  which  are  negligible.  Note  that  NB  can 
be  either  even  or  odd  for  the  magnetic  dipole  case  because  of  the 
symmetry  of  its  excitation  about  z  =  0,  and  NB  must  be  even  for  the 
electric  dipole  case  because  of  the  anti-symmetry  of  its  excitation. 


Fig.  A.l.  The  current  distribution  on  the  wire  induced  by  a  magnetic  dipole 

K.  u  ,  tor  r  =  0,001  and  d  =  0.1'.  (a)  The  real  part  of  total 

— k  R 

current . (h)  the  imas  inarv  part  ot  total  current .  (el  The  real  part 
of  evanescent  current.  (<1)  The  imaginary  of  evanescent  current. 
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Table  A.1.  Amplitude  of  outward  traveling  current  1^  on  the  wire 

induced  bv  a  magnetic  dipole  K.  ■  u  or  bv  an  electric 

— x 

dipole  I.  u  ,  for  r  =  0.001 >. 
v  B 
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APPENDIX  B 

COMPUTER  PROGRAM  LISTING  AND  DESCRIPTION 

The  computer  program  used  to  calculate  the  current  on  an 
inf initely-long  wire,  the  current  elements  of  equivalent  dipoles 
for  a  small  aperture  of  arbitrary  shape,  and  an  equivalent  network 
for  the  TEM  mode  on  the  wire,  is  described  and  listed  in  this 
appendix.  This  program  consists  of  subroutines  IMP,  CPI.,  SICI,  MUL, 

MUI.l,  GAUSS,  function  P  and  a  main  program.  A  subroutine  POL  used  to 
calculate  polarizabilities  of  a  small  circular  aperture  is  also  in¬ 
cluded. 

(A)  Subroutine  TM P 

Subroutine  IMP(NB,  K'L,  RB,  D,  Z)  computes  the  generalized  [Z] 
matrix  for  the  win'.  The  input  variables  are  defined  as 

NB  =  the  number  of  electric  current  expansion  functions  on 
the  wire . 

LI.  =  the  length  (in  wavelengths')  where  the  evanescent  current 
are  considered  to  exist. 

RB  =  the  radius  (in  wavelengths)  of  the  wire. 

D  =  the  distance  (in  wavelengths)  from  the  wire  to  the 
screen . 

The  output  is  stored  in  a  two-dimensional  array  Z.  Minimum  allocations 
are  given  by 

Complex  7.  (NB,  NB) 

DO  loops  10  and  20  compute  the  [7.1  elements  representing  the  coupling 
between  the  evanescent-current-segment  and  the  outward  traveling  currents. 


n  n 
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The  elements  for  the  coup]  Lngs  among  the  outward  traveling  currents 
are  then  computed.  This  subroutine  cal  Is  function  P  and  subroutine 
SICI. 


Listing  of  IMP 


c; 

C  IMF' {FIND  7  MATRIX:  J1=CEXP(.JK*Z)  .  J(NB)=CEXP(-JK*Z)  »J(2)~  J(NB-l) 

C  PULSE  EXPANSION  AND  IMPULSE  TESTING (EXCEPT  PULSE  TESTING  WHEN 

C  CALCULATE  Z(M.1>  AND  Z(M.NB) 

C  THE  •DELTA"  FUNCTION  OF  D(J1)/DZ  t  D(JNB)/DZ  (POINT  CHARGE 

C  IS  APPROXIMATED  BY  ‘PULSE"  CENTERED  AT  Z=0 

C  WL.:WIRE  LENGTH  FOR  HIGH  MODE  CURRENTS  ( J  ( 1  )  -  J(NB-1)> 

C  NB=  t  EXPANSIONS 

C  RB-RAD I US  OF  WTRF 

D-D 1ST .  BET.  WIRE  X  GND  PLANE 

SUBROUT  I  Nl  1  MP  <  NB  .  Wl  .  RB  .  D  .  7  ) 

COMF1  FX  7  <  NR  • NF< )  *  7MN  (  0  >  .P..IN.E 
DIMENSION  D7(f.> 

COMMON .'KKK/AK*  PI 
COMMON/  INK/  IN 
Al  B-  Wl.  /FL  OAT  (  ?*  (  NR  -?  )  1 
C 

C  ETND  Z(M.N>.  M.N-C?.  NB-  1  ) 

C  M-?  .  1  HE  1ST  SUBSECTION  ON  Wl  .  SO  M-l  IN  MO.N-l  IN  NO 

C  <MM.MO.MP>  «  ( NM  *  NO  *  NF‘  1  ARF  THE  POINTS  ON  SUBSECTION  M  %  N 

C 

DP-  D*? . 

NM  -  0 
NO-  1 
NP  ? 

DO  1  O  h-  '’.NR  -  1 
MOOJUM  I  i  1 
MP-MOn 
MM  MO  1 

P7  (  1  >  APB  (  f  I  OAT  (  MO  Nfl  >  #  Al  Bl 
D  7  (  ?  1  ADS  (FI  OAK  ME  *>'11*01  B) 

D7(  <  >  -  APS  <  1  I  OA  T  t  MP  NM  >  *At  B  > 

D 7  (  4  1  -ABS  (FI  OA  T  <  MM  Nf- 1  *AI  B  > 

D7  •  f.  >  -  ABS  (I  I  OAT  (  MM  WM  >  *AI  B  1 

do  :>o  ,ii.'. 

SO  7MN  '  I)  P  (  Al  B  .  D/<  ,n  .PR1  p(  A!  B.  D7  (  .1)  .  D?  > 

7<  M.  >  (  0  .  .  1  .  1  *740.  tf  T** ?*4*AI  B**?*7MN(  1  ) 

7(M,;m=:./im.'m  (0,  .1  ,  i*AO*(  7MN(0>  7MN(3>  -7MN ( 4 > +ZMN ( 5 >  ) 

7  (  ?  .  M  1  --  7  (  M  .  ?  ■> 

DO  TO  ,N  I.MH-1  M 
7(  m+.i.:m  .  I  >  -  7  <  M  .  ?  1 
7  (  ?+.l  .Ml  I  >  •  7  (M  •  .  1 
CONTINMI 


XO 
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o  n  n 
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FIND  Z  <  M  *  1 ) =7 (l.M).Z(M* NR )  =Z  (  NR  » M )  »  BY  USING  PULSE  TESTING 
OF  .1  ( M )  .  M=  ( 2  •  NP  - 1  ) 


40 

C 

C 

c 


A I  S=ALB 

022=112**2 

L  = ( NB- 1 )/2+l 
DO  40  H = 2 .  L 

70  - WL*FLOAT ( 2*M - NB- 1 > /Ft  OAT ( 2*N B -4 ) 

DZZ-ARS  <  ZO > 

ZNN <  1  )  240  .  *P J  *AI  P*  ( f ■  (  At.  B . DZZ » RB  > -P (  ALB »  PZZ » D2 )  > 
7P-ABS ( ZO+AL S > 

ZM=ABS ( ZO -AI...S  ) 

ZMN ( 2  >  (0.1 . )*AO.*(F(ALS»ZP.RB> -P  ( ALS  »  7P .  D?  > 

1.  -P  (  AL  S .  ZM .  RB  >  +P  (  Al. S .  ZM  •  D2 )  > 

Z(M.l  )  ZMN < 1 ) FZMN ( 2  ) 

7  (  M .  NB  )  --7MN  <  1  )  -ZMN(2) 

Z(l»M)-7(M.t > 

Z  (  NB  •  M  )  -7.  (  H » NB  > 

M  I  NB-  M+l 

/(Ml  ,  1  ) -2 ( M • NP ) 

7  (Ml  .NP) -  7 ( M • t > 

7  (  J  .Ml  >  =  Z ( Ml » 1  > 

7  (  NP  «M1  )  -~  7  (  M 1. « N  B  ) 

CONTINUE 

FIND  7 (  I • 1 ) ~7 ( NB . NB ) •  Z (I  . NB >  -Z ( NB  .  1  ) 

CALI  SICKS!  .Cl  .RP*2*PT  > 

CAI  I  .  S  T  C  T  <  S 2  .  C 2  .  n 2 * 2 *  P  T  ) 

ZMN  (  1  >  =  <  o  .  .  I  .  >  *A0 . * ( P ( AL  S . 0 . . RB  >  -P  ( Al  S  .  0  .  >  B2  >  ) 

ZMN ( 2  >  =60 .*(- C 1  PC 2  I (  O  .  .  1  .  >  *  (  S 1  -S2  )  ) 

7(1.1)  -  ZMN i  I  ) 

7(1  .  NP  '<  -  7MN(  1  -i  I  ZMN (2  1 

7  (  NR  »NB )  ■■•  /  (  1  .  I  > 

7  (  N  P  •  1  )  -  7  (  1  .  N  P  ) 

RETURN 
I  NB 
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where  X  is  the  input,  nr.  SI,  (  !  .ire  tin  outwits.  This  subroutine  is 
described  in  [ 1 h  ]  . 

Complex  function  P  t  \!  ,  /. ,  \l.)  evaluates  the  scalar  Green's  func¬ 
tion  .  defined  in  (t-a).  ii.e  i  v.i!  uit  ion  of  this  function  is  described 
in  | 1 0  j  .  The  input  variables  ar. 

Al.  =  i  ,  the  half  n.-th  ( in  wavelengths)  of  integration  for  .. 

B 

7.  -  the  lift  erence  (in  wave  1  enut  hs)  of  the  z-coortl  inatc  be¬ 
tween  tb.e  field  i'  ini  and  the  midpoint  of  source  element. 

XT  =  ,  the  tr.ansversi  eo.-rd  inate  (in  wavelengths)  of  the 

fie!  d  po  i  •’  f  . 


The  output  is  !' . 
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Listing  of  ST Cl 


subroutine:  sici<si>cifX) 

Z=AB3 ( X ) 

IF  ( Z~4 • ) 1  tit  A 
1  Y=(4.-Z)*(4.+Z) 

SI=X*< ( <  <  <1 ,753141E-9*Y+1 . 5A8988E-7 > *Y+1 .374168E-5) 

1  *Y+A.939889E~4)*Y+1. 9A4882E-2  >  *Y+4 . 395509E-1 ) 

CI  =  (  ( 5 . 7721 5AE-1  +AL.0G < Z )  >/Z-Z*<  (  (  <  ( 1 . 38A985E-10*Y 

1  +1 .58499AE -8>*Y+1 . 725752E-A > *Y+1 . 185999E-4 > *Y+ 

2  4 . 990920E-3 )  *  Y+l .315308E-1  > >*Z 
RETURN 

4  SI=SIN(Z> 

Y-COS  <  Z ) 

Z-;4  .  /Z 

U=  <<<(<(  < ( 4 . 048069E--3*Z-2 . 279143E-2 >  *Z+5 .  515070E-2 ) 

1  *Z~7. 261 642E~2>*Z+4, 987716E-2 > *Z~3 . 33251 9E-3)*Z- 

2  2 . 31461 7E-2  >  *Z~1 . 1 34958E-5 ) *Z+A . 2500 1 1 E-2 ) *Z+ 

3  2 . 5  8  3  9  8  9  E  - 1 0 

V=< ((((<( < ( -5 . 108699E- 3*Z+2 .8191 79E-2 ) *Z~6.537283E~2) 

1  *Z  +  7 . 902034E-  2  )  *Z--4 . 40041&E-2  )  *Z-7 . 94555&E-3  )  *Z+ 

2  2 . 60 1 293E  --2  >  *Z~3 . 764000E-4  >  *Z-3 . 1 224 18E-2  >  *Z- 

3  6 . A46441E- 7  >  *Z  +  2 . 500000E-1 
C I  --Z*<SI*V-  Y*U> 

SI=-Z*<SI#U+Y*V>+1 .570796 

RETURN 

END 


Listing  of  P 


C 


C 


COMPLEX  FUNCTION  P(AI...Z,ZL> 

COMPLEX  PO » ,JK 
REAL  II  .  12  *  1 3  « 1 4  *  LR «  KL 
COMMON/KKK/AK » P I 
COMMON/ JKK/JK 
RESORT <ZL**2+Z*#2> 

PO=CEXP <  - JK*R ) / ( 4 . #P1 ) 

I  r  <  R  .IT.  ( 1 0*AL. )  )  GO  TO  30 

FOR  R  >*  10* AL  r  FIND  P=*P<  A0>  A1 »  A2  »  A3»  A4  > 

KL-AN#AL. 

I  R-AI  /R 

ZR2-  <  Z  /  R  )  *  *  2 

ZR4  -  (Z/R)**4 

A00---1  .+3.  *ZR2 

A01  3 . -  30 . *ZR2+35 . *ZR4 

AO  1  .  +LR**2/A .  *A00+l.R**4/40  .  *A01 

A1--LR/6.*AOO  +  LR**3/40.*A01 

A2--ZR2/A . -LR**2/40 .*(1.-12 . *ZR2+15 . *ZR4 ) 

A3  --LR/60  .  *  <  3  .  *ZR2-5  .  *ZR4  ) 

A4=ZR4/120. 

P-PO/R* <  A0+  (  0  »  *  1  .  >*KL*Al  +  KL.**2*A2+(0.  1 1  .  ) *KL**3*A3+KL**4*A4  ) 
GO  TO  40 


.J 
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c 

C  FOR  R  10*ALf  FIND  P=P < 1 1 f I2f 13 f 14 )  I 

C 

30  Z 1 -Z+AL 

2.2  -7.  -AL 

G 1  -  S  Q  R  T  (  Z I  *  *  2 + Z 1  *  *  2 ) 

G2 -  SORT <  ZL  **2+Z2**2> 

C 

c  for  z  al: 

c 

1 1  -ALOG (  (Zl+Gl  >*<  -  Z2+G2)/Zl.**2> 

I F  (  /  .  LE  .  AL.  >  GO  TO  50 
C 

C  FOR  Z  >  AL : 

c 

1 1  -ALOG ( (Zl+Gl >/<Z2+G2> ) 

50  12-2 . *AL 

I3=Z1/2.*G1  -Z2/2.*G2+ZL**2/2.*I1 
I4<l.*AL*ZL**2+<2.*AI**3+4.*AL*Z**2)/3. 
P-P0/(2.#AL.>*<I1-JK*(  I2-R#I1 ) 

1  -AK**2/2.*( I3~2.*R*I2+R#*2*I1 ) 

2  +(0. f 1 . >*AK**3/A.*< 14-3. *R*I 3+3. *R**2*I2-R*#3*I 1 ) ) 

40  RETURN 

END 


(C)  Subroutines  ('!’].  .ini;  Pi'1 


Subroutine  CPI  (NR,  K'l.,  D,  T)  computes  the  coupling  matrix  [T] 
of  a  wire  lo  a  small  aperture  whose  magnetic  current  M  is  defined 
as  a  two-term  expansion.  The  inputs  NB,  Wl,  and  D  are  defined  in  the 
description  of  TMP.  The  outputs  are  stored  in  a  2  y  NB  matrix  T. 

The  minimum  allocation  is  given  hy 


Complex  T  (d,  N'B) 

Subroutine  POP  (PA,  AIM,  AI.E)  computes  the  magnetic  polariza¬ 
bility  and  electric  polarizability  for  a  small  circular  aper¬ 

ture  of  radius  r,.  The  input  is  RA  =  r,,  and  the  outputs  are  ALM  =  a 

A  A  mx 


and  AI,F.  = 


e 


n 
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Listing  of  CPI. 


c: 

r  cpi  :  find  the  t  matrix  for  a  small,  arperture 

C  DJDIST.  BETWEEN  WIRE  %  ARPERTURE 

C  WLtWTRE  LENGTH  FOR  HIGHER  ORDER  MODE  CURRENT 

C  THE  CENTER  OF  ARPERTURE  IS  AT  (O.O.O) 

c  nbj*  expansions: 

c 

SUBROUTINE  CPL  (NB.WL  .  D.T) 

COMP!  EX  T ( 2  *  NB ) » CRO » CR 1 . CR2 » CR 

COMPLEX  IK 

COMMON/KKK/AK  . PI 

COMMON/.  IKK/ JK 

AL.B=WL/F  I  OAK  2*  <  NB - 2  )  ) 

1)2  ~D**  2 

PD=1  ./(2.*PI*D> 

PI 2-PI *2. 

CRO=CEXP < - JK*D> 

T  (  1 .  1  >=-CRO*PD 
T  (  7 . 1  )  -CRO/  ( P 1 2*1*2 ) 

7  (1. NB  ) -T (1.1) 

T  <  2 . NR  >  ~-T ( 2 « 1 ) 

DO  10  N=2.  (NB-1.  )/2+l 
C 

THE  1ST  SUBSECTION  IS  N=2.  SO  2N-3  IN  ZO*. 

Z0=  FLOAT  (  2 * N - 3  )  *ALB-WL./2  . 

R--SORT  ( D? *70**2  > 

CR-CEXP  ( -•  JK*R ) 

T  ( 1  ,  N  )- ~D*ALB/P  I  *  <  JK/R**2+1  .  /R**3  )  *CR 
Rl-SQRT  ( 1*2+  (  701  AL  B  )  **2  ) 

R 2= SORT ( D 2  + ( ZO  -ALB) **2 ) 

C  R 1  "  ( 1  .  /  R  1**3  I  ..IK  7  R 1  *  *2 )  *CEX P  (  -  J K * R 1  ) 

CR2~  <  1  .  /R2**3+JK /R2**2  )  *CEXP  <  - ,JK*R2  ) 

T ( 2  » N  > - D /P 1 2* ( CR 1  CR2 ) 

N1=NB-N+1 
T  ( .1. .  N 1 )  -  T  ( 1 » N ) 

T  (  2 .  N 1 )  -  ~  T  (  2  •  N  ) 

10  CONTINUE 

RETURN 
END 


Listing  of  POI. 


SUBROUT  I NL  POI  *  PA  .  Al.  M  » ALE  ) 
ALM=4 , /3 . * PA * *3 
ALE” A I  M/2 • 

RETURN 

END 
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(D)  Subroutines  Ml'L,  Ml.'  1,1  and  CAV S S 

Subroutine  Mil.  (I.,  M,  X,  A,  H.  C)  computes  the  product  of  matrices 
A(L,M)  and  B(M,X>  and  stores  the  result  in  C’,(L,N).  The  inputs  are 
integers  L,  M,  X  and  matrices  A,  B.  The  output  is  matrix  C.  Minimum 
allocations  are  given  by 

Complex  A(I.,M>,  P.(M,X),  C(J  ,X) 

Subroutine  Mi'Ll  (I.,  M,  A,  B ,  C)  computes  the  product  of  matrix  A(L,M) 
and  vector  B(M),  and  stored  the  result  in  vector  C(L).  The  inputs  are 
integers  I.,  M,  matrix  A,  and  vector  B.  The  output  is  C .  Minimum  alloca¬ 
tions  are  given  by 

Complex  A  ( I . , M 1  ,  iUMl ,  (  (Ll 

Subroutine  CAl'SS  (X,  A,  B,  ITS,  ISVl  is  used  to  solve  a  linear  system 
equation  A  X  =  B  with  complex  coefficients  by  the  method  of  Gaussian 
elimination.  The  inputs  are  the  coefficient  matrix  A(N,X),  constant  vector 
B(X),  and  a  small  constant  FI’S.  The  output  I  Si-.'  is  1  if  the  absolute  value 
of  the  pivot  of  column  is  larger  than  ITS,  and  0  otherwise.  The  solution 
X  is  stored  in  B  as  an  output .  Minimum  allocations  are  given  by 
Complex  A(X,X1,  BfX* 


I 


Listing  of  MUL 
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SUBROUTINE  MUL. (L.M.N.A.B.O 
COMPLEX  A<L.M>.B<M.N)»C<L.N>»W 

DO  20  I  =  1.L 
DO  20  K  =  1  » N 
W=<0,0> 

DO  10  J=1»M 

10  W- “A(I» J)*B( J.K>+W 

20  C  <  I .  K  ) =W 

RETURN 
END 


Listing  of  Mi'Ll 


SUBROUTINE  MUL 1 < L » M . A . B » C > 
COMPLEX  A  <  L  .  M  >  .  » <  M  >  .  C  <  L.  > 

DO  10  1=1. L 
C  <  I )  =  (  0  ,  .  0 .  ) 

DO  10  J=1»M 

10  C(I)=A(I» J)*B< J)+C<I> 

RETURN 
END 


I 


j 


o  ©  o  © 
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Listing  of  GAUSS 


SUBROUTINE  GAUSS < N> A  * B , EPS » I SW > 
COMPLEX  A<N»N),B(N),C,T 
NM1 =N~ 1 
DO  10  K=1 , NM1 
C~ <0,0, 0.0) 

DO  2  I ~K f  N 

IP ( CABS (  A  <  I  >  K  ) ) . LE . CABS <  C ) )  GO  TO  2 
C~A  < I  *  K  > 

1 0  - 1 

2  CONTINUE 

I F ( CABS <  C ) . GE .EPS  >  GO  TO  3 

ISU-0 

RETURN 

3  IF(IO.EQ.K)  GO  TO  6 
DO  A  J~K,N 

T-  A  <  K  r  ,J  ) 

A  <  K  »  J  )  •-  A  ( 1 0  *  J  ) 

A  A  (  1 0  »  J  )  “  T 

T"-B  ( K ) 

B ( K ) (10) 

B ( 10) ~ T 
6  K  P.1  ••-KM 

C-l./C 
B(K)=B(K>*C 
no  10  J  -  K  F'  1  •  N 
A  <  K »  J  > =A ( K  *  J ) *C 
HO  20  I=KF'l  f  N 

A  <  .1  >  J  )  -A  ( 1  ,  J  >  -A  <  I ,  K  )  *A  (  K  t  J ) 

B(  J)-B(  .J)-  A(  J»K>*B(K) 
B(N)--:B(N)/A<N»N) 
no  40  K-l , NM1 
I  -  N  -  K 

o  <o,o,  o.o) 
i  p  i  ~ 1  *  i 
no  50  J-TPt  ,N 
C  '  C+A ( I »  J ) *B ( J ) 

8  < I ) =B< I )  C 
ISW=1 
RETURN 
END 
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(E)  Main  Program 

The  main  program  computes  the  coefficients  I  which  determine  the 

electric  current  on  the  wire,  coefficients  V  which  determine  the  cur- 

n 

rent  elements  of  the  equivalent  dipoles  for  the  small  aperture,  and 
the  equivalent  circuit  elements  for  the  transmission  line  mode  on  the 
wire.  It  calls  subroutines  IMF1,  POL,  GPL,  MUL,  MIJL1  and  GAUSS 
which  are  described  previously. 

The  input  data  arc  assigned  to  Lhe  program  according  to  the  DATA 
statement,  an  example  is 

DATA  N'B/1 2  / ,  WL/0.1/,  RB/0.001/,  RA/0.01/,  D/0.01/ 

where  NB,  U'L,  RB,  RA  and  I)  are  defined  previously.  The  solutions  are 

stored  in  an  output  file  named  'OUT'  whose  I/O  unit  is  assigned  to  be 

21  according  to  the  OPEN  statement. 

In  DO  loop  30,  current  distributions  on  the  wire  and  the  aperture 

are  computed  twice.  The  first  time  is  for  TF.M  wave  incident  from 

region  b,  and  the  second  time  for  plane  wave  incident  from  region  a. 

o  u 

During  calculating  circuit  elements  7^  and  Z^,  a  test  is  made  to  avoid 
negative  value  of  ReCz'j)  or  Re(Z^).  If  the  negative  value  is  negligible, 
a  value  of  zero  is  assigned  to  it.  If  it  is  a  large  negative  number,  an 
error  message  "***  STOP,  NEED  LARGER  WL  6.  NB  ***"  appears  in  the  output 
and  the  program  stops.  Usually,  NB  and  U'L  must  be  chosen  in  such  a  way 
that  WL/(2NB-4)  is  limited  in  fhe  range  between  RB  and  D,  and  ICL/D  2  10* 

In  addition,  NB  must  be  an  even  integer  number,  and  the  larger  P  is,  the 
larger  WL/D  is. 


n  n  r>  n 
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Minimum  alloc.it  Lons  are  given  by 


Complex  Z(NR,  NB)  ,  \Z(NA,  NA)  ,  i  T  ( N'A ,  Nil). 

CZ(NB,  NB),  DGZ(NB,  NB)  ,  VI  (NB)  , 
VM(NA),  ZZ (NB ,  NB),  VI (NA) , 

T(NA,  NB) ,  TT(NB,  NA) 


Listing  of  Main  Program 


COMPLEX  Z<12*12>*YZ(2»2)»YT<2»12)»GZ<12rl2)» 

1  DGZ(12*12)»VI(12) 

COMPLEX  UM(2) »7Z< 12 » 12) rYI<2> »AI1(2) » A 12(2) fZl t Z2»YTI (2) 

COMPLEX  F5»V1.U2 

COMPLEX  JK «  T  <  2 »  I  2  >  ?  T  T  ( 1 2  •  2  > 

EQUIVALENCE  (GZ.ZZ) 

COMMON/KKK/AK.PT 
COMMON./ JKK/JK 

DATA  NB/1  2/  -  WL../0 ,  1  /  »  RB/O  .  001  /  » RA/O  .  01  /  » D/0 .01/ 

DATA  ALAMPA/1 .0/ 

C 

C  NE<--#  OF  EXPANSI0NS<2  TEM f  NB-2  PULSE  EXPANSIONS) 

C  Ui: LENGTH  WHERE  HIGHFR  ORDER  MODE ( NON -TEM )  CONSIDERED 

C  RB  ’.RADIUS  OF  WIRE 

RAJ RAD I US  OF  CIRCULAR  APERTURE ( USED  FOR  SUB.  POL) 

D. DIST »  BETWEEN  WIRE  AND  APRETURE 
OIJTJ OUTPUT  FILE 

OPEN  (  UNIT  -21.  *EILE=  '  OUT  '  ) 

PI  =  3. 1415927 
AK-2 . *PI/ALAMDA 
JK  -  ( 0 . . 1 , )*AK 
NA-2 

1  CONTINUE 

WRITE  (21 ? 100)  NA  t NB ?  RA » RB  *  WL  t D 

C 

C  FIND  IMPEDANCE  MATRIX  OF  WIRE * Z ( NB > NB ) i 

C  INVERSE  OF  ADMITTANCE  MATRIX  *  YZ(2»2)5 

C  COUPLING  MATRICES  T  <  2  »  N.B  >  *  TT  <  NB  *  2  ) 

C 


CALL  I MP ( NB  *  WL  r RB  »  D  »  Z ) 
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C  OBTAIN  It  1  FOL  E  POLARIZABILITIES  ALM. ALE 

C  HERE  WE  CONSIDER  A  CIRCULAR  APERTURE  OF  RADIUS  RA 

CALL  POL.  < RA  » ALM  .  Al E ) 

El~8.*PI/3 
E2-120 . *PI 
E3=240.*FM**2*ALM 
E4=ALE*AO. 

YZ  <1*1 ) =F 1 /E2- <  0 . , 1 . )/E3 
YZ<  1  . 1  )~1 .  /  Y  7  <  1.  .  1. ) 

YZ<2.2)=PI**2/45. f (0. .1 . >/E4 
Y  Z  ( 2 . 2 )  1  .  /  Y  Z  (  2 »  2  > 

YZ  ( 1  »  2  )  --  (  0  .  .  0 .  ) 

YZ ( 2 • 1 ) = < 0 . . 0 . ) 

CALL  CPI.  <NB.UI  .D»T> 

DO  10  M--1.NB 
DO  10  N  1 » 2 

10  TT  <  M .  N  >  -  -  T  <  N  •  M ) 

C 

c  matrix  CALCULATIONS  I 

c 

CALL  MUI .  f  NA  *  NA . NR , YZ  » T . YT  > 

CALL  MUI.  <NFi.NAfND.TTf  YT.GZ) 

DO  20  M=}.NFs 
DO  20  N -■  1  .  NP 

20  DGZ ( M • N  ) -GZ <  M «  N  )  -  7  <  M .  N ) 

C 

K“.l  ?  TF  M  INCIDENT  <1+  =U!  K-2J  PLANE  WAVE  INCIDENT  (E  =1  ) 
C 

70  ■■60  ,  ♦AI..OG  <  2  .  #D/RB ) 

El-  F'T*D 
WRITE <  21 .1101 

DO  30  K-1.7 
VM<  1  )-j.  ,/E1 

VM  <  2  >  =•  ~  <  0  *  .1  ,  )*70/<<50.*El ) 

I E ( K  .EQ.  n  GO  TO  35 
WRITE (21 .120) 

Vh( 1 )=2. /< I  20. *PT ) 

VM  (  2  >  ~  <  0  ♦ . -2  * >/A0. 

35  CALI  MUI.  1  <NA.NA.YZ.VM.YI ) 

CALI  MUI.1  f  NB.NA.TT.YI  »VI ) 

DO  40  M-l.ND 
DO  40  N--1.ND 
40  ZZ ( M . N  >  =DGZ  <  M . N  > 

CALL  GAUSS ( ND . ZZ . VI . IE- 1 1 . ISW) 

IE  <  ISW  •  F.G .  1)  GO  TO  45 
TYPE  101 

101  F  ORMAT  < '  I SW-0 . STOP ' ) 

STOP 

45  WRITE ( 21 • 130 ) 

WRITE  (2.1  . 140)  <VI<M)  .M-l  »NFO 
All <K> -VI < 1 > 

AI 2  <  K  >  “VI  (  NF< ) 

CALL  MUL1 (NA.ND.YT, VI .YTI > 

DO  <50  M-1.NA 

60  V M<M) ■ YT (M)-YTI <M> 

WRITE <21 . 150) 

VM  <  2  >  =  <  0 . .  - 1 . ) #VM <  2 ) / 60 . 

WRITE <21 .200) <  VM<  M ) . M-l .NA) 


30  CONTINUE 
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C 

C  FIND  EQUIVALENT  NETWORK  ?Z1*Z2*V1*V2  (NORMALIZED  TO  ZO  ) 

C 

E5=2 . +AI 1 <1 >+A12( 1 > 

Z1=-(AI1  (.1  >+AI2<  1 )  )/E5 

Z2=2 .  *  ( 1  . +AI2  < 1 1 ) /E5/(  All < 1 )-AI2< 1 ) ) 

C 

C  TO  AVOID  ANY  POSSIBLE  SMALL  ERROR! SMALL  NEGATIVE  RESISTANT) 

C  DUE  TO  THE  VERY  SMALL  ROUNDOFF  ERROR  IN  REAL ( I ( 1 ) > , REAL ( I ( NB > > 

C 

RZ  -REAL (Zl  > /CABS ( Zl  ) 

IF (  RZ  .GT.  0.)  GO  TO  25 

IF<  ABS(RZ)  .GT.  l.E-02  )  GO  TO  27 

Z1=(0. » 1 . >*AIMAG(Z1  > 

25  RZ=REAL  <  Z2  )  ./CABS  (  Z2  ) 

IF  (  RZ  .GT.  0.  ) GO  TO  37 

IF ( ABS ( RZ )  .GT.  l.E-02)  GO  TO  27 

22- ( 0 , *1  , ) *A I MAG ( 72 ) 

37  CONTINUE 

E5=.t . fZl LZ2 

Vl=- E5*A!1 (2>+Z2*AI2(2) 

V2«F5*A 12  <  2 ) -72* A 11(2) 

WRITE (21  *1 70 )  71  *  72  *  VI  *  V2 

Zl..  =  ZO 

WRITE (21.  *190)  ZL 

WRITE (21. » 1  60)AI.1 .  (?)  *AT2(2) 

GO  TO  1.7 

27  WRITF(21 *210) 

1.00  FORMAT !//,2X*  '  NA-  '  *  14*  1  X*  '  NB-  '  *14*  IX*  '  RA=  '  *F6.4*1X* 

1  '  RB=  ' *  F6 . 4  *  1 X  *  '  Wt.  =  '  *  FA  » 4  *  1 X  *  '  D=  '  *  FA .  4  > 

110  FORMAT (/ *2X* 'TFM  WAVE  INCIDENT') 

120  FORMAT!/, 2X* 'PLANE  WAVE  T NCI DENT *E=1 ' ) 

130  FORMAT (2X* 'El  EC.  CURRENT  ON  WIRE* 1ST  X  LAST  ARE  TEM  AT  Z=0~  ,0+') 

140  FORMAT (2X*2E) 

150  FORMAT! /*2X* 'EQUIVALENT  DIPOLE  KL *  IL  =') 

160  FORMAT !/,2X* 'TEM  CURRENT  AT  Z=0-*0+  J ' , / * 2X , ' I 1 = ' * 2E * 2X * ' 12= ' * 2E ) 

170  FORMAT ( /, 2X  * ' EQUIVALENT  NETWORK ( NORMAL IZED  TO  ZO)J  ' , / * 2X * ' Zl  =  ' , 

1  2E*2X*  '72=' *  2E  * / * 2X  * 'Vl=' ,2E*2X* 'V2=' *2E) 

190  FORMAT ( / • 2X  * ' LOAD  IMPEDANCES  AT  TWO  PORTS  ZL  =  ZO=  '*F> 

200  FORMAT ( 2( 2X* 2E ) ) 

210  FORMAT ( 2X. '**#  STOP*  NEED  LARGER  WL  %  NB  *** ' ) 

17  CLOSE  <  UN T  T="2 1  •  F I L  E= '  OUT  '  ) 

STOP 

END 
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